#ifndef MOL_HPP
#define MOL_HPP

#include <iostream>
#include <math.h>
#include <vector>
#include <fstream>
#include "Eigen/Dense"

using namespace std;
using namespace Eigen;

double inftynorm(VectorXd v)
{
    double m = 0;
    for (int i = 0; i < v.size(); i++)
    {
        if (m <= abs(v(i)))
            m = abs(v(i));
    }
    return m;
}
#endif